Adsorption and Depletion of Polyelectrolytes from Charged Surfaces 
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Mean-field theory and scaling arguments are presented to model polyelectrolyte adsorption from 
semi-dilute solutions onto charged surfaces. Using numerical solutions of the mean-field equations, 
we show that adsorption exists only for highly charged polyelectrolytes in low salt solutions. Simple 
scaling laws for the width of the adsorbed layer and the amount of adsorbed polyelectrolyte are 
obtained. In other situations the polyelectrolyte chains will deplete from the surface. For fixed 
surface potential conditions, the salt concentration at the adsorption-depletion crossover scales as 
the product of the charged fraction of the polyelectrolyte / and the surface potential, while for a 
fixed surface charge density, cr, it scales as a^^^ f'^^^ , in agreement with single-chain results. 

PACS numbers: 82.35.Gh, 82.35.Rs, 61.41. -|-e 
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I. INTRODUCTION 

The phenomenon of adsorption of charged polymer 
chains (polyelectrolytes) to surfaces has generated a great 
deal of interest due to its numerous industrial applica- 
tions and relevance to biological systems. The theoretical 
treatment is not yet well established because of the mul- 
titude of length scales involved, arising from different in- 
teractions: electrostatic interactions between monomers 
and counter-ions, excluded volume interactions and en- 
tropic considerations. Furthermore, when salt is added 
to the solution, the interplay between polyelectrolytes 
(PEs) and salt ions as well as the ion entropy has to be 
taken into account. 

The adsorption of PE chains onto charged surfaces 
has been addressed theoretically in several models in 
the past. They include among others: solutions of lin- 
earized mean-field equations Q, Q, 0. S IE B , n umer- 
ical solutions of full mean-field equations |3, M, llOl- var- 
ious scaling theories for single-chain adsorption [ill , 
and formulation of a phenomenological criterion describ- 
ing the adsorption-depletion transition from charged sur- 
faces 

Q El H Other app roaches employed multi- 
Stern layer models where a discrete lattice 
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is used and each lattice site can be occupied by either a 
monomer, a solvent or a small ion. The electrostatic po- 
tential can then be calculated self-consistently together 
with the concentrations of the monomers and counteri- 
ons. 

In this article we re-examine the mean-field equations 
describing the PE adsorption and their numerical solu- 
tions, with specific emphasis on the adsorption-depletion 
transition. The present paper can be regarded as an ex- 
tension of Ref. Q. It agrees with the previously ob- 
tained low-salt adsorption regime but proposes a differ- 
ent interpretation of the high-salt regime. We find that 
the high-salt adsorption regime of Ref. is pre-empted 
by an adsorption-depletion transition, in analogy with 
single-chain results. The mean-field equations and their 
numerical solutions are formulated in Sec. ^1 some sim- 
ple scaling relationships in Sec. IIIII and the adsorption- 
depletion transition in Sec. IIVI A general discussion and 
comparison with other models are presented in Sec. 



II. THE MEAN FIELD EQUATIONS AND 
THEIR NUMERICAL SOLUTION 

Consider an aqueous solution of infinitely long PEs, 
together with their counterions and an added amount of 
salt. Throughout this paper we assume that both the salt 
ions and counterions are monovalent. Let 0(r) = \/c(r) 
be the square root of c(r), the local monomer concentra- 
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tion, a the monomer size and / the charge fraction on 
each PE chain. Also let (pb = ^/cb be the square root of 
the bulk monomer concentration Cb, and ipir) the elec- 
trostatic potential. The mean-field free energy can be 
obtained either from phcnomcnological or field theoreti- 
cal approaches: 



/ dr(/poi + /k 



/el) 



(1) 



/poi = k^T 



i«(04-04b)-M(0^-0g))(2) 



(c' - cl)] (3) 



/ci = (c+ -C- + fcj,') - — |VVr 



(4) 

While the full details can be found in Refs. |i|9|,ll0j, here 
we just briefly explain each of the terms. The first term of 
/poi accounts for chain elasticity, the second describes the 
excluded volume interaction between monomers, where 
V is the second virial coefficient. The third accounts 
for the coupling with a reservoir with bulk polymer con- 
centration (j)^ = Cb and chemical potential /i. The /ion 
contribution to the free energy takes into account the 
entropy of small ions and their chemical potential /x^. 
Lastly, /el is the electrostatic free energy. Its first term 
is the interaction energy between the electrostatic poten- 
tial and the charged objects; namely, the small ions and 
monomers. The last term is the self-energy of the electric 
field -^/drlVV'l^ 

Minimizing the free energy with respect to ^, (j), c'^ , , 
and using the bulk boundary conditions: ijj {x oo) = 0, 
oo) = 0b, = Csait + /0b s-iid c+ = Csait, the 
profile equations of Ref. jSj] are reproduced: 



C = (Csalt + /0b) ' 



(5) 
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s\nhpei> + (0^/0^^"^ - /02) (7) 
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— V20 = «(03-020) +/?/eV;0 



(8) 



where j3 = l/k^T is the inverse of the thermal energy 
ksT. Equations Q and © show that the small ions 



obey Boltzmann statistics, while Eq. ([7|l is the Poisson 
equation where the salt ions, counterions and monomers 
can be regarded as the sources of the electrostatic poten- 
tial. Equation (|SJ) is the mean-field (Edwards) equation 
for the polymer order parameter 0(r), taking into ac- 
count the excluded volume interaction and external elec- 
trostatic potential ip{r). 

The adsorption onto a flat, homogeneous and charged 
surface placed at x = depends only on the distance x 
from the surface. In this case the above equations can be 
reduced to two coupled ordinary differential equations. 
Defining dimensionless variables 77 = 0/0b and y = f^eip, 
Eqs. iQ and © then read: 



d^ 
dx^ 



K sinh y + kl^ (e^ 



6 da;^ 



vfl^l {v^ -v) + fyv 



(9) 



(10) 



where = (STr^eCsait)"^^'^ is the Debye-Hiickel screen- 
ing length, determining the exponential decay of the 
potential due to the added salt. Similarly, k^^^ = 
(47rZB0b/) " determines the exponential decay due 
to the counterions. The Bjerrum length is defined as 
Ib — e'^/eksT. For water with dielectric constant e — 80, 
at room temperature, Ib is equal to about 7A. Note that 
the actual decay of the electrostatic potential is deter- 
mined by a combination of salt, counterions, and polymer 
screening effects. 

The solution of Eqs. ^ and (|10|) requires four bound- 
ary conditions. Two of them are the boundary values in 
the bulk, X 00: 77 (x — > 00) = 1 and y {x —* 00) = 0, 
while the other two are the boundary conditions on the 
x — surface. In this article we use either constant 
surface charge density (Neumann boundary conditions) 
or constant surface potential (Dirichlet boundary con- 
ditions). For the former, dy/dx\x=o = —'ii^ue/ ekBT — 
— 47r(cr/e)ZB, where a is the surface charge density. For 
the latter, the surface potential is held fixed with a value: 
J/ (0) = j/s. The other boundary condition for the poly- 
mer concentration is taken as a non-adsorbing sur- 
face. Namely, (0) = 0. Note that far from the surface, 
X ^ 00, Eqs. and IjlOfl already satisfy the boundary 
condition: y = Q and (j) = fph (or 77 = 1). 

Equations and H10|) are two coupled non-linear dif- 
ferential equations that do not have a known analytical 
solution. The numerical solutions of these equations for 
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low salt conditions were presented in Ref. 8] and are re- 
produced here on Figs. 1 and 2, using a different numer- 
ical scheme. The numerical results have been obtained 
using the relaxation method j3| based on a lineariza- 
tion procedure done on a discrete one-dimensional grid. 
Then, the equations are transformed to a set of algebraic 
equations for each grid point. The sum of the absolute 
difference between RHS and LHS over all grid points is 
minimized iteratively until convergence of the numerical 
procedure is achieved. 

In calculating the numerical profiles of Figs. 1 and 2 
we assume positively charged polymers and a constant 
negative surface potential. In Fig. la the reduced elec- 
trical potential y = Peij) profile is shown as a func- 
tion of the distance from the x = Q surface. Similarly, 
in Fig. 2a the monomer rescaled concentration profile 
c(a;)/cb, is shown. In both figures a constant surface 
potential boundary condition is imposed. The differ- 
ent curves correspond to different surface potentials y^, 
monomer charge fractions / and monomer size a. From 
the numerical profiles of the electrostatic potential and 
monomer concentration it can be clearly seen that there 
is a distinct peak in both profiles. Although they do 
not occur exactly at the same distance from the surface, 
the corresponding peaks in Fig. 1 and 2 vary in a similar 
fashion with system parameters. The peak in the concen- 
tration (Fig. 2) marks a PE accumulation at the surface 
and is regarded as a signature of adsorption. The peak 
in the potential (Fig. 1) marks an over-compensation of 
surface charges. At the peak of i^ix), the electric field 
vanishes, E = —dip/dx = 0, meaning that the integrated 
charge density from the surface up to this distance ex- 
actly balances the surface charge. 

III. SCALING ESTIMATE OF THE 
ADSORPTION LAYER: COUNTERION ONLY 
CASE 

So far numerical solutions within mean-field theory, 
Eqs. lini and (|10|l . have been described. We proceed 
by presenting simplified scaling arguments, which are in 
agreement with the numerical mean-field results. Note 
that the treatment here does not capture any correla- 
tion effect which goes beyond mean-field. The concept 
of polymer "blobs" can be useful in order to describe PE 
adsorption, where such polymer blob can be regarded as 
a macro-ion adsorbing on a charge surface. The blob 



size is determined by taking into account the polymer 
connectivity and entropy as well as the interaction with 
the charged surface. A single layer of adsorbing blobs 
is assumed instead of the full continuous PE profile as 
obtained from the mean-field equations. Therefore, the 
blob size characterizes the adsorption layer thickness. 



A. Fixed Surface Charge Density 

The two largest contributions to the PE adsorption free 
energy are the electrostatic attraction with the surface 
and the chain entropy loss due to blob formation. 

For simplicity, the electrostatic attraction of the 
monomers with the surface is assumed to be larger than 
the monomer excluded volume and monomer- monomer 
electrostatic repulsion. With this assumption the chain 
has a Gaussian behavior inside each surface blob, D ~ 
ag^l"^ ^ where g is the number of monomers in a blob of 
size I?, as is shown schematically in Fig. 3. The entropy 
loss of the chain balances the surface-monomer attrac- 
tion. As a result the blob attraction with the surface is 
of order fcsT. It is now easy to get an estimate of the 
blob size D: 

fg-e^D~kBT (11) 

e 

g~ (?B/aa)"'/3 , (13) 

using a rescaled surface density a = \a/e\. These results 
are in agreement with those describing the statistics of 
single-chain adsorption [REt|. 

The assumption that the electrostatic attraction to 
the surface is larger than the monomer-monomer electro- 
static repulsion and excluded volume can now be checked 
self-consistently, yielding two conditions: a ^ /a^^ and 
/» j;3/(ai°aZB). 

The average monomer concentration (per unit volume) 
in the adsorption layer c^, is the blob concentration in 
the adsorption layer no, times the number of monomers 
per blob g, yielding €,„ = n^g- It is now possible to get 
an estimate of the blob concentration per unit volume in 
the adsorption layer, no, by assuming that the adsorbed 
layer neutralizes the surface charges up to a numerical 
prefactor of order unity. Hence, no — a/Dfg. This 
assumption, which is in agreement with our numerical 
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solutions, leads to 
no 



Ba (14) 
Cm = nog^ilBa-'r'ay/'^a'^'r'/' . (15) 



Equation p4|l is just the Graham equation |21| relating 
the surface charge density with the counterion density 
at the surface vicinity. The only difference is that the 
counterions are replaced by the charged polymer blobs. 
Furthermore, Eq. ifT^ is in accord with the results of 
Ref Q. 

The total amount of PEs in the adsorption layer is F ~ 
CmD = a/ f. In other words, the overall polymer charge 
in the adsorption layer (up to a numerical prefactor) is 
/F ~ a. This is just another way to phrase the charge 
neutralization by the PEs mentioned above. 



B. Fixed Surface Potential 

Using the boundary condition of a fixed surface poten- 
tial V' = V's , the scaling laws for D and g can be obtained 
in a similar fashion as was done in Sec. IIII Al Alterna- 
tively, one can (in the absence of salt) relate the surface 
potential to the surface charge density by ips — crD/e. 

The adsorption energy of a blob of charge gfe onto 
a surface held at potential "03 is just gfeips- Requiring 
that this energy is of order of fcsT we obtain in analogy 
to Eqs. (HU-linj): 



D ~ ag'^ ~ 



a 



vfWs 
1 

~ flVs 



(16) 
(17) 

(18) 



Together with the neutralization condition uq 
a/iDfg) it yields: 



no 



nog ' 



(19) 
(20) 



Note that the above results are in accord with the ones 
previously derived in Ref. . 

Just like in Sec. III. A the self-consistent check can be 
repeated here for the dominance of the surface-monomer 
interactions, yielding \ys\ 3> f^^^l^'^a"'^^^ and / ^ 
This condition has been verified, in addition, 
by examining numerically the mean-field adsorbing pro- 
files. 



The overall charge of the polymer in the adsorbed layer 
is then: 



F - 



\y. 



3/2 jl/2i 









dip 




IbD 


dx 


x=0 



-1/2 



(21) 
a (22) 



which again verifies that the adsorbed amount scales like 
the surface charge. 

The numerical results of the mean-field equations for 
constant surface potential condition and in the low salt 
regime (cgait — 0.1 mM) are consistent with this scaling 
picture, as can be seen in Figs. 1 and 2. In Fig. lb the 
rescaled potential y/\ys\ is plotted in terms of a rescaled 
distance: x/D, with D taken from Eq. H17|l . In Fig. 2b 
the concentration profile is rescaled by Cm , Eq. (|20|l , and 
plotted in terms of the same rescaled distance x/D. The 
figures show clearly data collapse of the two profiles, in- 
dicating that the characteristic adsorption length D is 
indeed given by the scaling predictions. Note that the 
agreement with the scaling argument occurs as long as 
the system stays in the low salt limit. The other limit of 
high salt is discussed next. 



IV. THE ADSORPTION - DEPLETION 
TRANSITION IN PRESENCE OF ADDED SALT 

The same numerical procedure outlined in Sec. II, 
is used to find when the chains stop adsorbing and in- 
stead will deplete from the surface. This is not a sharp 
transition but rather a crossover which is seen by cal- 
culating numerically the PE surface excess, as depicted 
in Fig. 4. The profiles were obtained by solving numeri- 
cally the differential equations for several values of / near 
the adsorption-depletion transition using a fixed surface 
potential boundary condition. For salt concentration of 
about c*a^n ~ 0.16 \ys \ J /{Ibo^) (solid line in Fig. 4), the 
figure show the disappearance of the concentration peak. 
Namely, a depletion-adsorption crossover. 

The dependence of F = Jg°° dx {(jP' — (pl) on Cgait and 
/ for constant surface potential is presented in Fig. 5. 
The place where F = indicates an adsorption-depletion 
transition, separating positive F in the adsorption regime 
from negative ones in the depletion regime. In Fig. 5a the 
dependence of F on / is shown for several salt concentra- 
tions ranging from low- to high-salt conditions. For low 
enough /, F < indicates depletion. As / increases, a 
crossover to the adsorption region, F > 0, is seen. In the 
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adsorption region, a peak in r(/) signals the maximum 
adsorption amount at constant Cgait- As / increases fur- 
ther, beyond the peak, F decreases as 1/ \fj . Looking at 
the variation of F with salt, as Csait increases, the peak 
in F(/) decreases and shifts to higher values of /. For 
very large amount of salt, e.g., Csait — 0.5 M, the peak 
occurs in the limit / ^ 1. In Fig. 5b, we plot F(csait) 
for several / values. The adsorption regime crosses over 
to depletion quite sharply as Csait increases, signaling the 
adsorption-depiction transition. The salt concentration 
at the transition, Cg^^j^, increases with the charge fraction 
/. The dependence of F on Csait and / for constant sur- 
face charge density is plotted in Fig. 6. Both salt and / 
dependences show a similar behaviour to those shown in 
Fig. 5 for constant surface potential. 

The numerical phase diagrams supporting the 
adsorption-depletion transition are presented in Fig. 7 
for constant surface charge conditions. The phase di- 
agrams were obtained by solving numerically the mean- 
field equations. We scanned the (/, Csait) parameter plane 
for 50 values of / between 0.01 < / < 1 (Fig. 7a) and 
the ((T, Csait) plane for 50 values of a ~ |c/e| between 
10-^ [A-2] < < 10--* [A-2] (Fig. 7b). From the log- 
log plots it can be seen that the adsorption-depletion 
transition is described extremely well by a line of slope 
2/3 in both Fig. 7a and 7b. Namely, at the transition 



''salt 



j2/3 fixed a and 



'salt 



^2/3 Jqj, fixed /. 



To complete the picture, the adsorption-depletion 
transition is also presented in Fig. 8 for constant sur- 
face potential. The phase diagrams are obtained by solv- 
ing numerically the differential equations. We scanned 
the (/, Csait) parameter plane for 50 values of / between 
0.01 < / < 1 (Fig. 8a) and the (ly^l ,Csait) plane for 50 
values of \ys\ between 0.1 < \ys\ < 1.0 (Fig. 8b). From 
the figure it is apparent that the adsorption-depletion 
transition line fits quite well a line of slope 1.0 in both Fig. 
8a and 8b plotted on a log-log scale. Namely, c*ait ~ / 
for fixed ys, and c*g^j^ ~ ys for fixed /. 

These scaling forms of c*^^ at the transition can be ex- 
plained using the simplified scaling arguments introduced 
in Sec. in. 



charges up to distances D from the charged surface. For 
high ionic strength solutions, the electrostatic potential 
at distance x for a charged surface can be approximated 
by the linearized Debye-Hiickel potential: 



y {x) = 'iiralBK 



(23) 



This is valid as long as the potential is low enough, y < 
1. The adsorption picture requires that the exponential 
decay of the potential will not vary substantially inside a 
region of size D comparable to the size of surface blobs, 
y{D) ~ j/s- Then, the exponential decay in Eq. I|23|) 
yields 



kD <1 



(24) 



Namely, the Debye-Hiickel screening length is smaller 
than the adsorption layer thickness, D. Using Eq. H12() 
this yields: 



Csait < o- 



~2/3^2/3;-l/3^-4/3 



(25) 



The crossover between adsorption and depletion will oc- 
cur when c*jju ~ (a^ f^l^^a^^y^"^ , in accord with Refs. 
0, 0, 01 , and with the numerical results discussed above 
and presented in Fig. 7. 



B. Scaling for Fixed Surface Potential 

For the boundary condition, ip ~ ipg, the potential 
decay from the surface can be approximated to be: 

y (x) = y.e-'^^ (26) 

and the same consideration as in Eq. H24|) and (|25|1 gives: 

\ys\f 



Csait < 



(27) 



Namely, we expect an adsorption-depletion transition to 



occur for 



'salt 



lUsl //(^Bfl^), in the case of a fixed sur- 



face potential. This supports the numerical results as 
presented in Fig. 8. 



V. DISCUSSION 



A. Scaling for Fixed Surface Charge 

If the blobs are taken as charged spheres, the mere 
existence of an adsorption process requires that the at- 
traction of the monomers to the surface persists for all 



We have presented numerical calculations of the mean- 
field equations describing the adsorption of PE chains 
onto charged surfaces, including multi-chain interactions. 
The main finding is the existence of an adsorption- 
depletion transition in presence of added salt or weakly 



charged chains. The numerical results are discussed in trie, 
terms of simple scaling arguments describing the adsorp- 
tion of PEs. The salt concentration at the adsorption- 
depletion transition scales like c^^^ ^ f \ys\ for fixed sur- 
face potential and c*^^^ ^ (/o')^^'^ for fixed surface charge 
density. Within the scaling picture, the condition for de- 
pletion is the same as for a single chain, in agreement 
with our mean-field solutions. 

We briefly summarize the main approximations of our 
mean-field and scaling results. A non-adsorbing surface 
is used as the polymer boundary condition. However, if 
the surface has a strong non-electrostatic affinity for the 
PE chains, the electrostatic contribution does not have to 
be the dominant one. The method also assumes Gaussian 
blobs within mean-field theory. In a more refined theory, 
excluded volume interactions as well as lateral correlation 
in the blob-blob interactions will alter the adsorption be- 
havior. When the surface charge (or potential) is high 
enough, the blob size D can become comparable with 
the monomer size a, and the PE chains will lay flat on 
the surface. Further investigations might be necessary to 
address in more detail the above points. It will also be 
interesting to extend our results to geometries other than 
the planar charged surface. 

Several authors have addressed the problem of adsorp- 
tion onto surfaces either of a single chain or mul- 
tiple chains |3| using similar arguments of blobs. In 
another approach, a Flory-like free energy Q was in- 
troduced using the assumption of a single characteris- 
tic length scale. The latter gave adsorption-layer scaling 
laws as in Eqs. (|17|) and (|20|l . but did not find the deple- 
tion criterion. Instead, an adsorption length scale and a 
characteristic concentration were predicted for the high- 
salt regime. We show here, using both numerical calcu- 
lations and scaling arguments, that the high-salt regime 
does not exist because it is preempted by a PE depletion. 
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Figure Captions 

Fig. 1 (a) Numerical profiles of the rescaled electrostatic 
potential y = fieip as function of the distance from 
the surface x using Eqs. and IjlUI) . and constant 
surface potential. The solid line is for a = sA, 
f — l,ys = —1.0, the dotted line for a — 5A, / — 
l,ys — —0.5, the dashed line is for a — lOA, / — 
1, j/s = —0.5, and the dash-dot line for a = 5A, / = 
0.1, j/s — —0.5. All profiles have Csait = 0.1 niM, 
= 10-6 A-3, V = 50A3, e = 80, T = 300 K. The 
profiles reproduce those of Ref. j^l using a different 
numerical scheme, (b) Same profiles as in part (a) 
but in rescaled variables: x/D and y/ |t/s|. 

Fig. 2 (a) The concentration profile c(a;)/cf, = (j)^{x)/(j)^ 
for the numerical calculations specified in Fig. 1. 
The profiles reproduce those of Ref. . (b) Same 
as in part (a) but in rescaled variables: x/D and 
c{x)/cm- 

Fig. 3 A schematic drawing of polyelectrolyte adsorption 
onto flat surfaces and formation of Gaussian surface 
blobs each of size D and having g monomers. The 
monomer size is a. 

Fig. 4 Numerical polyelectrolyte concentration profiles 
exhibiting the transition from adsorption to deple- 
tion. The dashed line corresponds to f — 0.12, the 
dot-dash line to / = 0.1, the solid line to / = 0.09, 
and the dotted line to / = 0.08. All profiles have 
\y^\ = 0.5, <j)l = IQ-^A-^, V = 50A3, a = 5A, 
Csait — 70 mM. The adsorption-depletion transition 
is found to occur for / = 0.09, corresponding to 
c:^,,^0.16\ys\f/lBa\ 

Fig. 5 (a) Surface excess of PE adsorption, T, as func- 
tion of the chain charged fraction /, for several salt 
concentrations: l.OmM (solid line), lOmM (dashed 
line), 0.1 M (dash-dot line), 0.5 M (dots), and for 
constant surface potential. As the salt concentra- 
tion increases, the peak in F shifts to higher / val- 
ues and disappears for Cgait = 0.5 M. The depletion- 
adsorption transition occurs for F = 0. (b) Surface 
excess as function of salt concentration, Csait, for 
several / values: f=0.03 (dots), 0.1 (dashes), 0.3 
(dot-dash), 1.0 (solid line). F is almost indepen- 
dent of Csait for low salt concentrations in the ad- 
sorption region. It is then followed by a steep de- 
scent into a depletion region at a threshold value. 



Other parameters used are: ys — —l-O, v — 50A^, 
= 10-6A-3, a = 5A, T = 300K and e = 80. 

Fig. 6 (a) Surface excess of PE adsorption, F, as function 
of the chain charged fraction /, for several salt con- 
centrations: 4.0 mM (solid line), 8.0 mM (dashed 
line), 21 mM (dash-dot line), 63 mM (dots), and for 
constant surface charge density. As the salt con- 
centration increases, the peak in F shifts to higher 
/ values and disappears for Csait = 63mM. The 
depletion- adsorption transition occurs for F = 0. 
(b) Surface excess as function of salt concentra- 
tion, Csait, for several / values: f=0.1 (dots), 0.2 
(dashes), 0.45 (dot-dash), 1.0 (solid line). F is 
almost independent of Csait for low salt concen- 
trations in the adsorption region. It is then fol- 
lowed by a steep descent into a depletion region 
at a threshold value. Other parameters used are: 
a/e = -10-4A-2, V = 50A3, ^2 ^ 10-6A-3, 
a = 5A, T = 300K and e = 80. 

Fig. 7 Numerically calculated adsorption-depletion 
crossover diagram for constant surface charge 
condition. In (a) the (/, Csait) parameter plane on 
a log-log scale while a = \a/e\ is held constant 
at CT = lO-^A-^. The full squares represent the 
lowest salt concentration for which depletion is 
detected. The least-mean-square fit to the data 
points gives a straight line with slope of 0.69 ±0.02. 
The figure shows that the numerical results agree 
with a 2/3 power law as predicted in Sec. IV. A, 
'^sait ^ /^^'^- In (b) the crossover diagram in cal- 
culated numerically in the (|(T/e|, Csait) parameter 
plane on a log-log scale, while / is fixed to be 
/ = 0.1 The least-mean-square line has a slope 
of 0.71 ± 0.02, showing that the numerical results 
agree with a 2/3 power law as predicted in Sec. 
IV.A, cU ^ a2/3. 

Fig. 8 Numerically calculated crossover diagram on a log- 
log scale for constant surface potential conditions. 
Notations and symbols are the same as in Fig. 6. 
In (a) the (/, Csait) parameter plane is presented 
for constant t/s = —1.0. The least-mean-square fit 
has a slope of 1.00 ± 0.02, in excellent agreement 
with the scaling arguments, Cg^it ~ /■ In (b) the 
(|ys|,Csait) parameter plane is presented, for con- 
stant / — 0.1. The least- mean-square fit has a 
slope of 1.04 ± 0.02, in agreement with scaling ar- 
guments, c*^n - |ys|- 
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